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ABSTRACT 

We identify a correlation space involving optical and UV emission line parameters as well 
as the soft X-ray spectral index that provides optimal discrimination between all principal 
classes of active galactic nuclei. Most of the sources in our three high quality data samples 
show a strong intercorrelation with narrow line Seyfert 1 (NLSyl) galaxies and steep spectrum 
radio galaxies occupying opposite extrema in the space. NLSyl sources show a clear continuity 
with broader line sources indicating that they are not a disjoint class of AGN as is sometimes 
suggested. We interpret the principal intercorrelation in the parameter space as being driven 
by AGN luminosity to black hole mass ratio (L/M oc Eddington ratio). Source orientation no 
doubt also plays an important role but it is not yet clear whether FWHM H/3 or CivA1549 
line shift is the better indicator. We tentatively identify two RQ populations: an almost pure 
RQ population A with FWHM<4000 and population B which occupies the same parameter 
domain as the flat spectrum RL sources. A possible interpretation sees Pop. A/NLSyl as lower 
mass/high accretion rate sources and Pop. B/RL sources the opposite. 



1. Introduction 

The search for correlations among observational parameters that describe AGN has been developing 
rapidly in the past few years. At the same time there has been a theoretical countercurrent that views 
all quasar spectra as remarkably similar. This erroneous view is due in part to the low resolution and 
s/n of much available spectral data. Blurred quasar spectra do look remarkably similar but data with 
s/n>20 in the continuum and resolution ~5A show striking differences from which a pattern is beginning 
to emerge. Some of the most promising correlations arose when high s/n optical spectra were published 
for the lower redshift (z ^ 0.5) sources in the Bright Quasar Survey (Boroson & Green 1992: BG92). A 
principal component analysis of the BG92 correlation matrix showed "eigenvector 1" correlations involving 
the widths and strengths of [OHIJA5007, as well as broad H/3 and FEli opt emission lines. More recently 
measures of X-ray luminosity, particularly the soft X-ray photon index, have emerged as a related part of 
the eigenvector 1 correlations (Wang, Brinkmann & Bergeron 1996). 

We report a study of correlations involving the best available data samples for low redshift (z ^ 1.0) 
Active Galactic Nuclei (AGN). We find a correlation space where we are able to discriminate between all 
of the major forms of AGN phenomenology (e.g. broad and narrow Line Seyfert 1 galaxies, regular and 

1 Department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487; giacomo@merlot.astr.ua.edu 
2 Department of Physics, University of Ljubljana, Jadranska 19, 1000 Ljubljana, Slovenia; tomaz.zwitter@uni-lj.si 
3 Osservatorio Astronomico di Padova, vicolo delFOsservatorio 5, 1-35122 Padova, Italy; marziani@pd.astro.it 
4 Instituto de Astronomia, UNAM, Mexico, DF 04510, Mexico; deborah@astroscu.unam.mx 



-2 - 



broad absorption line (BAL) quasars, steep and flat spectrum broad line radio galaxies). We refer to it as 
the Eigenvector 1 (hereafter El) space reflecting its partial roots in the BG92 study. The three principal 
"orthogonal" correlates involve measures of: 1) low ionization broad line width - full width at half maximum 
of H/3 (FWHM H/3bc); 2) ratio of line strengths - the continuum normalized ratio of the optical FEii opt 
and broad H/3 emission line strengths (R,Feii= W(Feii A4570 blend)/W(H/?Bc) and 3) X-ray continuum 
strength - the soft X-ray photon index (r so f t ). Parts of this correlation space have been discussed for almost 
ten years (e.g. Boroson & Green 1992 (PG92); Boiler et al. 1996; Marziani et al. 1996 (MS96); Wang et al., 
1996; Laor et al. 1997) but the pieces have not previously been united in this way. We present the basic El 
phenomenology in this paper. In a following paper we will consider the principal physical drivers of the El 
correlations. 



2. The El Correlation Space 

Figures [l] a,b,c show the 2D projections of the El correlation space. AGN plotted in the figures come 
from three data sources: 1) the BG92 sample of 87 mostly radio-quiet (RQ) AGN (16 radio-loud: RL), 2) 
the MS96 mixed sample of 21 RQ and 31 RL sources as well as 3) a sample of 24 new sources with matching 
HST UV archival and optical ground-based spectra (8 RL). The data sources combine the most complete 
sample of AGN with high quality spectral data (BG92) and two overlapping samples of comparable quality 
data with matching optical H/3 and UV CivA1549. Soft X-ray photon indices are available for a large 
fraction of the above samples (Wang et al. 1996; Brinkmann et al. 1997; Siebert et al. 1998; Yuan et 
al. 1998). Values determined with hydrogen column density N# as a free parameter were preferred when 
available. The hard X-ray photon index provides an alternative, but less sensitive, El correlate (Piccinotti 
et al. 1982; Brandt et al. 1997). Our total sample includes 128 sources (45 RL) with optical measures, of 
which, 76 (39 RL) have UV measures. Mean error bars (2a) are indicated with: i) average errors for sources 
in the middle of the diagrams for Fell A4570 and H/3 measures as well as ii) a median value for all r so f t 
measures. See BG92 and MS96 for detailed discussion of reduction and analysis procedures. 

RQ sources are indicated by solid symbols in the figures while open and crossed circles indicate flat 
(core-dominated) and steep (lobe-dominated) spectrum RL (as defined by Kellermann et al 1989) sources 
respectively. The majority of sources show a well defined correlation in Figures |l] a,b,c. One can make a 
first order interpretation of El correlations in terms of RL vs. RQ differences. Sources that fall in the large 
FWHM H/3 and large R Fe ii region (like PG 0043+039; Turnshek et al 1994) are apparently very rare and 
pathological (possibly associated to mixed Starburst/ AGN properties). Fig. [l] 1 a,b,c indicate that these 
two populations show a clear separation with RQ having, for example, a mean FWHM(H/3bc) ~ 2300 
km s _1 less than RL sources (a few RL sources with FWHM(H/3bc) ~ 11-20000 km s _1 are not shown in 
the figures). An alternative interpretation of El sees two populations of RQ AGN: 1) population A (filled 
boxes) which shows little overlap with the RL domain (65% of the BG92 RQ sample) and 2) population 
B (filled circles) that occupies essentially the same El domain as the RL sources ( 25% of the BG92 RQ 
sources). Fig. Q 1 d shows a plot of W(Feii A4570) vs. FWHM(H/3 BC ) for the BG92 sample that provided 
the motivation for the population A-B concept. Rather than a correlation this plot shows two disjoint 
populations of sources with no correlation within either population. Correlations appear when RfoII is 
used instead of W(FEii opt ) or W(H/3bc)- We introduce the population A and B distinction for the sake of 
argument but suggest that it may have more fundamental significance. Table 1 presents mean parameter 
values (and associated sample standard deviations) in arbitrary bins of FWHM(H/3bc) that allow one to 
compare RL vs. RQ as well as Pop. A vs. Pop. B sources. 
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It is not easy at this time to determine what role sample selection biases play in El but it is clear that 
BG92 selection techniques favored RQ sources with narrow Balmer line profiles. Recent work indicates that 
selection on the basis of soft X-ray excess favors them even more strongly (e.g. Moran et al. 1996; Grupe 
et al. 1999). At the same time core-dominated/flat spectrum RL sources are no doubt over-represented if 
their emission is beamed. This study makes no claim to completeness. Beyond BG92 our sample selection is 
driven by the availability of i) high s/n and resolution ground based data and ii) matching HST UV spectra. 
Table 2 presents the results of a correlation analysis for the three principal El and related parameters 
(the Pearson's correlation coefficient rp is reported along with the probability P of a chance correlation in 
a sample of N objects). As the figures and Table 2 also suggest, the strongest El correlations are found 
for the RQ sources. The highest ranked El correlation (r so f t vs. FWHM H/?bc) is present only in the 
RQ sample which is not surprising since RL sources show no evidence for a soft X-ray excess. The nearly 
identical rp for the ALL RQ and Pop. A RQ samples suggest that the principal correlations are driven 
by Pop. A. RQ Pop. B and RL sources show the same El parameter domain and neither sample shows 
evidence for significant internal correlations. The latter result is due, at least in part, to: a) the absence of 
a soft X-ray excess in these sources and b) the weakness of FEii opt emission in these relatively broad-lined 
sources. If the minimum detectable FEii opt emission for sources with profile widths FWHM(Feii A4570) w 
FWHM(H/3 B c) «1000 km s _1 is W(Feii A4570) w 10-15 A then we estimate that it will be « 25 and 40 A 
for FWHM(H^bc) « 5000 and 10000 km s" 1 , respectively. The message here is that values of R Fe ii ~ 0.2 
are very uncertain even with the best available data. 

There are several lines of evidence that suggest continuum luminosity is uncorrelated with El parameter 
space. 1) Optical luminosity appears in the second orthogonal eigenvector identified in BG92, 2) BAL 
sources occupy a similar domain to the NLSyl but they are "displaced" in: a) X-ray luminosity - BALs 
are X-ray quiet and b) optical luminosity where they are on average about 4x more luminous than the 
NLSyl sources in BG92. 3) RL and Pop. B RQ sources occupy the same parameter domain but the RL 
sources show ~ 5 x higher optical continuum luminosities which may in part reflect the presence of beamed 
sources in the RL sample. 4) We also find no evidence for a correlation between El coordinates and radio 
continuum luminosity for RQ sources (see Kellermann et al. 1989; Falcke et al. 1996). This last result tells 
us that RL sources are a distinct AGN population that shows fundamental differences in BLR structure and 
kinematics. El indicates that the weak radio emission from RQ sources is unrelated to the RL phenomenon. 

El also shows evidence for a parameter space separation between steep and flat spectrum RL sources. 
This shows up in Table 2 as a possible RL correlation between FWHM H/3bc and RfcII- Steep spectrum RL 
sources represent the opposite extremum from NLSyl while flat spectrum sources are more similar to RQ 
sources. The tendency for steep spectrum RL sources to show the broadest Balmer profiles and the weakest 
FEii op t emission is confirmed in two large (with some source overlap between themselves and our RL 
sample) surveys (Brotherton 1996; Corbin 1997). Siebert et al (1998) provide evidence that this separation 
may also be present in the X-ray spectral index. NLSyl are a RQ and population A extremum with the 
narrowest Balmer line profiles, strong FEii opt emission and a strong soft X-ray excess. NLSyl show a clear 
continuity and correlation with broader line Seyfert 1 galaxies in all measured parameters. This challenges 
the idea that they represent a unique or disjoint AGN population. BAL quasars (4 RQ BAL in the BG92 
sample) occupy an El domain that is similar to the NLSyl. Other (BAL) studies (e.g. Boroson & Myers 
1992) have also suggested that BAL quasars show Balmer line FWHM ;$ 3000 km s _1 and moderate to 
strong Rpciinieasures. 

Much less clear is whether RQ population B sources represent a disjoint AGN population or show 
a smooth continuation of the population A correlations. Figure b for example shows a breakdown of 
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the r so f t - FWHM(H/3) correlation at the nominal boundary between population A and B. Whatever the 
relation between the two RQ populations, their line profiles show striking differences. Figure || shows a 
comparison of the H/3bc an d CivA1549bc l me profiles for prototype NLSyl source I Zw 1 and NGC 5548 
which is a typical broad line Seyfert galaxy. While differences in the Balmer profiles are striking, the most 
impressive difference is the apparent kinematic decoupling of the CivA1549 and H(3 profiles in I Zw 1. Our 
previous work (MS96) suggested that properly NLR corrected CivA1549 in RQ sources (see Sulentic & 
Marziani 1999) is always blueshifted relative to H/3 and the AGN rest frame. Fig. || shows the CivA1549bc 
line shift vs. FWHM(H/3 B c) diagram. We have normalized the CivA1549 shift by W(CivA1549) in this 
plot because we see a complementary trend in El for W(CivA1549) to be smallest in the NLSyl population. 
Figure || confirms that essentially all RQ sources show a CivA1549 blueshift while RL sources show equal 
red and blueshifts with amplitudes generally less than ±10 3 km s _1 . Correlations of ClvA1549 shift with 
RfcII and r so f t (see Table 2) indicate that it is likely to be an important El correlate. We will explore it 
further in succeeding papers. 



3. From Observational to Physical Parameters 

El correlates AGN spectroscopic data in a way that removes much of the apparent "randomness" of 
line properties. It also redefines input parameters for photoionization and kinematical models. We have only 
begun to explore the El parameter space and and the physics that drives it. Accumulating evidence from 
numerous studies suggests that the correlations in El involve at least two principal independent parameters: 
the AGN luminosity to black hole mass ratio and the source orientation. The role of Fe abundance, disk 
magnetic fields and black hole angular momentum are beyond the scope of this introduction to El. 

Available evidence suggests that the central source luminosity to black hole mass ratio L/M may be 
systematically higher in Pop. A sources (i.e. a higher accretion rate). A role for L/M is supported by: a) 
the existence of a soft X-ray excess, particularly in NLSyl, that has been related to a higher accretion rate 
(Pounds et al. 1995) and b) the thermal signature of the disk itself (Page et al. 1999; Puchnarewicz et al. 
1998). If the soft X-ray excess is the high energy end of a thermal signature of the accretion disk then a 
steep rising blue optical continuum (in low z sources) can be interpreted as the low energy end of the "big 
blue bump". This would explain why BG92, with quasars selected on the basis of U-B colors, would favor 
detection of NLSyl. CivA1549 shows the largest blueshifts in NLSyl sources, and the lowest CivA1549 W 
values, both consistent with the idea that the line emission arises in a disk wind (e.g. Murray & Chiang 
1997; Bottorff et al. 1997). The observations are consistent with the idea that L/M is only high enough 
in RQ sources (especially Pop. A) to trigger a radiation pressure driven outflow. The Balmer profiles in 
Pop. A sources are most easily interpreted as arising in an optically thick illuminated accretion disk (see 
e.g. Sulentic et al. 1998). RL sources show Balmer lines with both large red and blue profile shifts. 

RL sources show no evidence for a soft X-ray excess which may reflect a real absence of this component 
or that it is dominated by a relatively hard component related to X-ray emission related to the radio 
loudness. We favor the former interpretation because most RQ population B sources also show a soft X-ray 
deficit. The absence of a soft X-ray component in RL+RQ pop B sources coupled with the absence of the 
CIV blueshift are consistent with weak or absent disk structure. The broad lines may arise in a biconical 
structure in many or all of these sources (Marziani et al. 1993; Sulentic et al. 1995). Evidence for a bicone 
origin of the BLR includes: a) sources with double peaked broad lines (e.g. Sulentic et al. 1995; Eracleous 
& Halpern 1998), b) sources with single peaked broad lines showing both large red and blue displacements 
(e.g. Marziani et al. 1983; Halpern et al. 1998), c) sources with a transient BLR (e.g. Storchi-Bergmann et 
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al. 1993; Ulrich 2000) and spectropolarimetry of a) and b) sources (Corbett et al. 1998). 

The El parameter space is measuring aspects of both the geometry and kinematic of the broad line 
region in AGN. The strength of the correlations and reasonable orthogonality of the parameters suggest 
that a better diagnostic space for AGN is unlikely to be found. Even in the preliminary presentation, El 
allows us to resolve some AGN conundrums. 1) RL sources are found to be an AGN population with 
fundamental geometrical and kinematic difference from the RQ majority. 2) NLSyl are found to be an 
cxtremum of the RQ population rather than a pathological or disjoint AGN population. 3) If the width 
of the Balmer lines bear any signature of source orientation then some NLSyl (at least, I Zw 1) are seen 
at or near pole-on (face-on accretion disk) orientation, while some BAL QSOs are likely to be misaligned 
NLSyl. 4) Population B RQ quasars with Balmer lines broader than w 4000 km s _1 may represent a distinct 
class (RL pre/post cursors) from narrower population A RQ sources. 5) CivA1549 line shifts indicate the 
presence in Pop. A (and possibly Pop. B) RQ sources of an ubiquitous disk outflow/wind. 

PM acknowledges financial support from the Italian Ministry of University and Scientific and 
Technological Research (MURST) through grant Cofin 98-02-32 and from the IA-UNAM where part of this 
work was done. 
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Fig. 1. — Panels a,b,c: The three principal correlation planes of Eigenvector 1. Solid symbols are radio-quiet 
sources (boxes: Pop. A, circles: Pop. B). Open symbols are radio-loud sources (open: core dominated, 
crossed: lobe dominated). Panel d: the BG92 sample plotted in the W(Fcii A4570) vs. FWHM(H/3 B c) 
plane. The dotted line at 2000 km s _1 separates NLSyl and the rest of the sample, while the dot-dashed 
line at 4000 km s^ 1 sets the limit between Pop. A and Pop. B. Error bars indicate 2a uncertainties for a 
"typical", non-pathological data point at FWHM(H/3bc) ~ 4000 km s -1 , R,F e ii~0.5, and error on r so f t equal 
to the median of errors for which r so f t was available. 
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Radial Velocity (km s _1 ) 



Fig. 2. — H/3 (right panels; solid line: H/3 B c) and CIVA1549 (left panels) broad line profiles for NLSyl 
source IZwf ("extreme" Pop. A, upper half) and Pop. B Scyfcrt 1 source NGC5548 (lower half). Note that 
the abscissa origin of I Zw 1 has been set by Hi 21 cm measures. [Oiii]A5007 lines show an appreciable 
blueshift with respect to the Hi reference radial velocity as well as to the peak of the H(3 profile. Dotted 
lines: residuals after broad component (thick line subtraction). Dashed line in I Zw 1 H/3 panel: residual 
spectrum after FEii opt subtraction. 
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Fig. 3.— Relationships between FWHM(H/3 B c) and CivA1549 FWHM profile shift (blueshifts are negative) 
normalized by the rest-frame equivalent width of CIVA1549, W(CivA1549). Error bars are indicative as in 
Fig. |l]and refer to W(CivA1549) w 75 A and shift « 400 km s -1 . 



TABLE 1 
Properties of AGN Populations 



AGN Class 




FWHM range 


N 


W(rL3) 


W(FeiiA4570) 


RfcII 


W([Om]A5007) 


r SO ft 


N r 






[km s _1 ] 




[A] 


[A] 




[A] 






RQ Pop. A 




0-2000 


20 


72 ± 28 


57 ± 18 


0.85± 0.29 


15.5 ± 12.0 


3.0 ± 0.44 


18 


RQ Pop. A 




2-3000 


19 


103 ± 36 


65 ± 27 


0.74± 0.42 


27.0 ± 38.0 


2.8 ± 0.37 


18 


RQ Pop. A 




3-4000 


14 


121 ± 39 


55 ± 17 


0.50± 0.25 


24.0 ± 23.0 


2.45 ± 0.31 


11 


All RQ Pop. 


A 


<4000 


53 


96 ± 39 


59 ± 21.5 


0.72 ± 0.35 


22.0± 27.0 


2.82 ± 0.43 


47 


All RQ Pop. 


B 


>4000 


27 


104 ± 33 


39 ± 27 


0.40 ± 0.29 


27.0 ± 25.0 


2.33 ± 0.28 


23 


All RQ 






80 


99 ± 37.5 


52 ± 25 


0.61 ± 0.37 


23.5 ± 26.0 


2.66 ± 0.45 


70 


CD RL 






19 


79 ± 29.5 


28 ± 17.5 


0.36 ± 0.16 


20.0 ± 11.5 


2.06 ± 0.73 


19 


LD RL 






2G 


88 ± 39 


16.5 ± 14.0 


0.22 ± 0.18 


25.0 ± 17.0 


2.25 ± 0.33 


23 


RL 




< 4000 


8 


68.5 ± 36 


24 ± 13 


0.41 ± 0.20 


13.0 ± 8.0 


2.06 ± 0.30 


8 


RL 




> 4000 


37 


88 ± 34 


21 ± 17 


0.25 ± 0.17 


25.0 ± 16.0 


2.18 ± 0.60 


34 


All RL 






45 


84 ± 35 


21.5 ± 16.5 


0.28 ± 0.18 


23.0 ± 15.0 


2.16 ± 0.55 


42 



TABLE 2 

El Sample Pearson's Correlation Coefficients 



Correlates 




ALL 




ALL RL 




ALL RQ 




P 


N 


rp 


P 


N 


rp 


P 


N 


rp 


P 


N 


rp 


FWHM(H/3) vs. R Fo „ 


125 


-0.50 


< 10- 5 


45 


-0.37 


0.01 


80 


-0.44 


4 x 10~ 5 


47 


-0.' 


FWHM(H/3) vs. T 


112 


-0.41 


1 x 10~ 5 


42 


0.11 


0.49 


70 


-0.63 


< io- 5 


47 


-0., 


r so f t vs. Rfcii 


112 


0.40 


1 x 10~ 5 


42 


-0.21 


0.18 


70 


0.58 


< 10- 5 


47 


-0., 


FWHM(H/3) vs. v r (Civ) 


70 


0.145 


0.23 


39 


-0.15 


0.36 


31 


0.19 


0.30 


18 


-0.( 


Rfcii vs. v r (Crv) 


70 


-0.52 


1 x 10~ 5 


39 


0.09 


0.59 


31 


-0.62 


2 x 10~ 4 


18 


-0., 


rsoft vs. v r (Civ) 


G5 


-0.20 


0.10 


37 


0.09 


0.60 


28 


-0.40 


0.35 


16 


-0., 


W(H/3) vs. W(Civ) 


72 


0.31 


0.009 


39 


0.23 


0.16 


33 


0.46 


0.007 


19 


0.4 



